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The preparation of materials with large nonlinear responses usually 
requires involved synthetic processes. Thus, it is very advantageous for 
materials scientists to have a means of predicting nonlinear optical 
properties. For example, correlating the strengths of the nonlinear 
properties of a family of compounds with different substituents at various 
bonding sites may be of great utility in selecting good candidates for 
experimentation . 

The prediction of nonlinear optical properties has to be addressed 
first at the molecular level and then as bulk material. At the molecular 
level, prediction of these properties must consider both conditions of 
static and dynamic fields. For relatively large molecules, two types of 
calculations may be used: the sum-over-states (based on a perturbative 
expansion in terms of molecular eigenstates in the absence of an external 
field) and the finite-field approach (which diagonalizes a molecular 
Hamiltonian with a static-field dipolar perturbation) . In the first case, 
the magnitude of the properties is obtained for a specific dynamic field, 
and the static field value can be extrapolated to zero frequency. In the 
second case, the nonlinear terms are obtained for static conditions and a 
correction for dispersion should be added if they are to be compared with 
experiments involving electromagnetic radiation. 

The finite-field method was selected for this research, because this 
approach is better suited for larger molecules. A simplified version of 
the sum-over-states approach is under development to predict dispersion 
effects. Previously, dispersion corrections were extrapolated from 
experimental data on molecules similar to the ones investigated. Finally, 
a molecular mechanics approach is being considered to predict interatomic 
effects on nonlinear properties. These effects have been estimated 
previously from empirical data on molecules similar to the ones under 
study . 

a) Progress in static-field calculations 

During the summer of 1992, a core correction to the third-order 
polarizabilities was implemented. The method to calculated static 
nonlinear polarizabilities is a modification to the method developed by 
Dewar and Stewart, 1984 for calculating molecular linear polarizabilities. 
Dewar and Stewart include isolated point charges (or "sparkles") in the 
description of the molecule, with no associated electrons or atomic 
orbitals. An appropriate selection of these sparkles simulates the effect 
of a uniform electric field on a molecule, allowing for the corresponding 
changes in polarization to be calculated directly. The all valence- 
electron MNDO Hamiltonian (Dewar and Thiel, 1977) is employed to obtain 
molecular polarizations in the presence of various static fields. Then, 
polynomial regressions of polarizations versus fields data are used to 
calculate all tensor elements of the static second and third-order 
polarizabilities. Since MNDO is a valence-electron model, a correction for 
the atomic cores was required for third-order polarizabilities. In second- 
order polarizabilities, the effect of induced polarization in 
centrosymmetric structures cancels since they involve squares of the field 
strengths . 
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A group of substituted benzene molecules was selected to parameterize 
atomic core corrections for carbon, hydrogen, oxygen, and nitrogen. 
Experimental data for these molecules was available from third-harmonic 
generation (THG) measurements (Cheng et al., 1991). in order to minimize 
solvent effects, the molecules chosen for comparison were those measured as 
neat samples. Since the THG measurements were done at 0.649 eV radiation 
the experimental values were corrected for dispersion effects. The 

correction for dispersion of the experimental values was based on the two- 
state anharmonic model: 
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V th 1 StatiC third " order Polarizability, y © is the measured 
t ird-order polarizability, © 0 is the resonance frequency (©. = 6.9 eV for 
benzene), and 0) is the laser frequency (Zhao et al., 1989). The 

!^!f imental data were also corrected for bulk effects, as 'follows- Two 
EFISH measurements on liquid and vapor benzene were used to obtain a liquid 
o vapor correction. The liquid-to-vapor correction was assumed to be 
applicable for all substituted benzene molecules. 

A plot of number of atoms in the molecules versus the difference 
between the predicted third-order polarizabilities and the corrected 
experimental value showed a correlation of 0.90 (Figure 1, solid line) 

from g th! nOEmali3!ation fact °« the same non-hydrogen to hydrogen ratios 

1984) a Dlot :7*: CtXO :\° n linear P° larizabil ities (Dewar and Stewart, 
1984), a plot of H units versus the difference between predicted and 

experimental values showed a correlation of 0.91 (Figure 1, dotted line) 

Thus the significant correlations indicate that core corrections fo^ ' 
third order polarizabilities were required. By performing multiple 
regressions with various groups of data, it was found that the core 
corrections for carbon and hydrogen were highly coupled. In order to 

correction f P roblem ' the same carbon to hydrogen ratio from the core 
° r i0n ° f linear Payabilities was used. Thus the following atomic 

»- s ;- 25 v 10-61 “• - 61 

SI T,hi 0 1 e fc K g 0-386 x 10 SI, and nitrogen = 0.931 x 10 -61 
oi. Table 1 shows the results. 

d *r° nd and th ^ rd -° rder Polarizabilities have been predicted for 

tffnnllT ^ 0mp ° Und3 cont aining pyrryl or substituted amines and a variety 
of function 31 groups. This work has been recently submitted for Y 

publication (Cardelino et al., 1992). 

Second-order polarizabilities have been calculated for nitro- and 

f^m°iif >St:LtUted qU ^ n ° lines and isoquinolines. The work was initiated 
from discussions with Dr. Wayne Hyde, a JOVE participant. Organic 

synthesis of these materials is in progress at the Northwestern State 

fTl^ Y t L ° UiSiana - A P a P er on this work will shortly be suited 
for publication (Moore et al., 1992) 

Second-order polarizabilities have been calculated for substitute 

nitrobenzene. The work wee initiated from discussions with another JOVE 

Dr - R °” ald Synthesis of these material is under 

at the New Mexico Highland University. Y 
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Predicted value THG 


10-61 

C m4 v-3 

H 

3.3 

3.5 

CH3 

4.2 

4.1 

0CH3 

4.8 

4.3 

NH2 

4.9 

4.8 

N(CH3)2 

6.3 

7.2 

CN 

4.7 

3.8 

C (0) H 

4.6 

4.8 


b) Progress in dispersion effects 

The basis for predicting the dispersion effects of nonlinear 
polarizabilities is time -dependent perturbation theory. At low intensities 
the P and y terms are constant . In order to detect second and third 
harmonics, the molecules must be exposed to the intense radiation of 
lasers. Under these circumstances, the nonlinear polarizabilities become 
dependent on the frequency of the incident light. 

During the summer of 1992, a method to estimate transition dipole 
moments was developed. To approximate the required transition dipole 
moments, the zero-differential-overlap method was applied to ground and 
excited states. Unoccupied molecular orbital wavef unctions from ground 
state calculations in the absence of external fields were used to construct 
various excited states. The dipole moments were separated into atomic 
point charge and hybridization contributions. 

Since the goal is to design a method suitable for predicting 
dispersion effects on large molecules, the relative importance of the 
perturbative terms will be considered. Inclusion of indirect or direct 
summations, the number of terms in each summation, corrected expressions 
for the excited states and molecular orbital energies will be studied. 

c) Progress on bulk effects 

In our previous work, the effect of intermolecular interactions on 
the second and third-order polarizabilities have been estimated from 
experimental data fit to solvent molecular refraction or solvent dielectric 
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on a parent compound was 


constant. In the case of neat samples, the effect 
assumed valid for the substituted molecules. 

Calculations performed on dimethylamino nitrobenzene have proven to 
be very sensitive to conformational isomerization. It is of great 
importance to know what the structure of the molecule is in the crystalline 
state. This problem will be examined using a molecular mechanics approach. 
Estimations of heats of sublimation and fusion, as well as hardness of 
crystals will also be done. Knowing these properties will help determine 
the feasibility of growing crystals of a given substance and its utility in 
various applications. 

During the summer of 1992, the molecular mechanics program MM3 was 
modified to be run in a 486 personal computer. 
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